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The surface of a topological insulator hosts a very special form of a quasi-two dimensional metallic 
system when it is embedded in a topologically trivial medium like the vacuum. The electronic 
properties of this unusual 2D metal are distinct in many aspects from both the conventional two- 
dimensional electron gas systems in quantum well heterostructures as well as those of a single layer 
graphene. In this paper, we study one of these distinct features i.e., the response of the electronic 
spins to an applied magnetic field perpendicular to the surface. We find an unusual behaviour of 
the spin magnetization and susceptibility as a function of both the magnetic field and the chemical 
potential for a generic topological surface. We propose that this behavior could be studied by the 
recently developed experimental technique called /3-NMR which is highly sensitive to the surface 
electron spins. We explain how this technique could be used to probe for spontaneous magnetic 
ordering caused by magnetic dopants or interactions discussed in the recent literature. 



I. INTRODUCTION 

Topologically protected electronic states residing on 
the surface of topological insulators (TI]P^form a unique 
2D metal distinct from those so far realized in solid state 
systems. Similar to the low energy electronic states in 
a single layer graphen^, the robust 2D metal on an TI 
surface has conic branches touching (in the absence of 
the intrinsic gap) at high-symmetry points in the first 
Brillouin zone. This resemblance in the energy disper- 
sion is responsible for some common properties between 
these metallic systems, such as the square-root depen- 
dence of Landau Levels (LLs) to the applied magnetic 
field. However, there are important differences between 
these systems as a result of the topological nature of the 
surface states in TIs. 

The low energy Dirac-like surface bands of a topo- 
logical insulator arise due to the mismatch of the bulk 
topological invariants on the two sides of the surface re- 
gardless of how constituent atoms have been arranged on 
the surfa ce as long as the bulk remains in the topolog- 
ical phas#^. In the case of a 'strong' topological in- 
sulator (STI) there is an odd number of Dirac points at 
any surface and, importantly, they are not spin degener- 
ate. Furthermore, they exhibit a unique spin- momentum 
lockin^^ as a result of the strong spin-orbit interaction 
in the underlying STI which is essential for the formation 
of the topologically non-trivial insulator phase with time- 
reversal symmetry. These properties make electronic sur- 
face states in a STI very distinct from those of graphene 
with a pair of spin degenerate Dirac-like bands near two 
special points in the Brillouin zone which exist as a re- 
sult of the unique arrangement of carbon atoms in the 
honeycomb lattice. 

One of the interesting aspects in which these 2D metal- 
lic systems behave uniquely is their magnetic response. 
Here, motivated by the recent progres in the experimen- 
tal methods and the importance of the spin susceptibil- 
ity for our understanding of electronic systems, we study 
the spin response of electrons on the surface of a STI 
and show that it exhibits interesting features even in a 



simple non-interacting limit. We find that the character- 
istic spin-momentum locking of these electrons leads to a 
unique spin response that is distinct from spin-degenerate 
systems like graphene and 2DEG. Instead of the oscilla- 
tory behaviour of the susceptibility as a function of the 
chemical potential found in spin degenerate systems, in 
STIs our work predicts a plateau-type behaviour which 
arises from the strong correlation between spin and or- 
bital degrees of freedom. In addition we find that the 
existence of a special LL with full spin polarization leads 
to a jump in the magnetization as the chemical potential 
crosses the energy of this LL. When an intrinsic mag- 
netic ordering is present there is also a jump in the mag- 
netization as a function of the applied magnetic field. 
Our results can ther efore assist in detection of such in- 
trinsic magnetizatiorP^^ which is known to have pro- 
found consequences for the nature of the surface state; 
the magnetized surface state of a TI is predicted to be- 
come a quantum Hall liquid with half-integer quantized 
Hall condu ctiyit\ *^^ * -^^ ' and many unusual^^^^^ and poten- 
tially usefuP^^ physical properties. 

There are various ways of measuring the weak mag- 
netization produced by electrons on a metallic surface. 
The SQUID scanning magnetometry is a highly sensi- 
tive probe which can detect tiny magnetization on the 
surfaces. However, it is a challenge to use this device 
to probe magnetization when there is a large applied 
magnetic field which interferes with the superconducting 
part of the SQUID and causes noise. A variant of this 
method using a superconducting pickup coil has been 
developed in order to study deHaas-van Alphen oscilla- 
tions of the 2D metallic systems in a large perpendic- 
ular magnetic field. High-sensitity micro- mechanical 
cantilever magnetometr};^^ is another way of measuring 
electronic magnetization, however, this method measures 
the magnetization of the whole sample and it is difficult 
to isolate the contribution from the electrons on a single 
surface. 

The nuclear magnetic resonance (NMR) is another 
powerfull experimental technique which can be used 
to study electronic spin magnetic response in a bulk 
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metalP^ Through the so-cahed Knight shift in the nu- 
clear resonance peak, it is possible to probe the spin 
part of the magnetic susceptibility of the electrons as 
they interact with the resonating nuclei in their prox- 
imity. Unfortunately, this method, in its conventional 
form, fails to be useful in very thin films and 2D metallic 
systems due to the limitation in the number of available 
nuclei and the resulting weakness of the signal. Thanks 
to the progresses made by experimentalists in control- 
ling and implementing high energy beams of unstable 
ions, an experimental technique, the so-called /3-NMR, 
has been recently developed to overcome the above limi- 
tations. Briefly, in this exotic variety of NMR, unstable 
radioactive ions such as ^Li and ^^Be are implanted in 
the sample surface. The nuclear spin precession signal 
is then detected through the products of the beta decay 
of the radioactive nucleus. Since the ion implantation 
depth can be controlled by tuning the beam energy it is 
possible to acquire information about the behaviour of 
the electronic spins in very thin metallic films. ^^^^^ Ex- 
periments are currently underway to study the surface 
magnetic response of STI crystals Bi2Se3 and Bi2Te3 us- 
ing /3-NMR.29 

The paper is organized as follows. First we introduce a 
simple model known to describe the electrons on the sur- 
face of a STI in the presence of a perpendicular magnetic 
field. We review the exact solutions of its eigenvalue 
problem whic h has been studied previously in various 
contexts. ^^'^ In section III, we calculate the spin mag- 
netization and susceptibility assuming that the chemical 
potential lies inside the gap between positive and neg- 
ative energy eigenstates. We discuss the magnetic re- 
sponse of the surface as one tunes the magnetic field and 
the chemical potential for various values of the intrin- 
sic gap. In closing, we explain how a /3-NMR experiment 
might be able to detect these effects through Knight shift 
measurements. 




FIG. 1: (Color online) The surface spectrum of a strong topo- 
logical insulator in the absence of magnetic dopants (Ao = 0), 
(a) in the absence of the external magnetic field, (b) in the 
presence of an applied perpendicular magnetic field. 



coupling of the magnetic vector potential, hk ^ hk-\- eA 
(electron charge — e). The other is the coupling of the 
spins to the magnetic field, the Zeeman effect, expressed 
as 5Hz = —gsl~iB^~^^ • 5, where Qs is the effective elec- 
tron gyromagnetic constant. In the bulk Bi2Se3 crystal 
Qs ^ 30 (Ref.[33]) although much smaller values have been 
reported for electrons near the surface. In our calcula- 
tions below we use two representative values, = 8 and 
gs = 30, which yield qualitatively similar results with 
some interesting differences. 

In the continuum limit, the leading order Hamiltonian 
describing these surface states in the presence of the ap- 
plied magnetic field, B = Bqz^ in the Landau gauge 
A = —{Boy^ 0), can be written in the following form 



(2) 



II. THE MODEL 

Many of the interesting features of the surface states in 
a topological insulator can be captured, at least qualita- 
tively, using a simple non-interacting Dirac Hamiltonian 



(1) 



where = (c| 



) and is the fermionic creation 



operator of the spin up (down) states with wave vector k. 
Ao is the intrinsic gap in the surface spectrum which 
might be nonzero when the time-reversal symmetry is 
spontaneously broken due to magnetic ordering. The lat- 
ter can arise due to the presence of magnetic dopants with 
spin S in the proximi ty of the surface exchange-coupled 
to the electronic s pin d 1 -*^^ ^ or as a result of electron- 
electron interactioiP^. 

Turning on a perpendicular magnetic field adds two 
terms to the above Hamiltonian. One is the minimal 



where {y) is the creation operator for the spinor 
mode extended along the x direction and localized at 
y. 1-L{kx^y) for Bq > {) \s defined as 



n{kx,y) 



where Cc = VFy^2eh\Bo\ and 



A = Ao 



(3) 



(4) 



The term in A proportional to the magnetic field is the 
Zeeman contribution, a^^ is the one-dimensional har- 
monic oscillator bosonic operator defined as 



elBn 



(5) 



Note that varying kx shifts the position of the localized 
state produced by the application of on the vacuum 



3 



state along the y direction. In the Bq <{) case, 1-L{kx^y) 
can be obtained from the one given for 5o > in Eq. ^ 
by exchanging the off-diagonal elements and replacing kx 
by its time reversed counterpart —kx- 

The eigenstates for Bq > {) and n > are given h^^^ 

(cos {3n/2)ifri-i{y - kxh^)\ 
-ism{dn/2)(pn{y - kxh )/ 



(sin {Sn/2)(pn-i{y - kxlB'^)\ 
2 ' 
I COS {6n/2)(pn{y - kxh ) / 

while for n = we have 

^k.Ay) =( ) • (8) 

\(po{y - kxh'^)/ 

In the above cos Sn = A/e^ and (/^^ are the 
one-dimensional harmonic oscillator eigenstates, i.e., 
a^acpn = TiLpri' We remark that the n = eigenstate 
in Eq. (|8| is very special since it is fully spin polarized. 
This will have important consequences for the magnetic 
response discussed below. 

The eigenvalues associated with cj)^ ^{y) eigenstates 
are given by 

4 = ±Vne2 + A2, n > (9) 

and for the fully spin-polarized n = eigenstates 

So = -sgn(5o)A. (10) 

Note that the form of eigenstates for Bq < {) is differ- 
ent from that given in Eq.(|6]|8| since the Hamiltonian is 
different in that case. 
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FIG. 2: (Color online) Spin magnetization 5Ms = Ms(T, /x) - 
Ms (0,0) in units of (xo- Tesla) as a function of the chemical 
potential for a nonmagnetic surface (Ao = 0) at Bo = 3.0 
Tesla, Qs = 8 (top panel) and 30 (bottom panel). /cbT = 
0.1, 1.0,5.0 meV (red, green, blue). 

M-r. = ^ • ^ (13) 



III. SPIN SUSCEPTIBILITY AND 
MAGNETIZATION 

The electronic magnetic moment due to spin is propor- 
tional to the spin operator 

fie = -7eS, (7e = -QsI^b/K) (H) 

Therefore, to calculate the spin part of the magnetic mo- 
ment for the eigenstates given in the previous section we 
only need to find the expectation value of the spin opera- 
tor. A straightforward evaluation using Eqs. ([6][8| shows 
that all the electronic states within the same LL con- 
tribute equally to the magnetization. For each of them 
we have 



Ml^ = Ml„ = (12) 



M.,o = -'-^sgn{Bo) (14) 

The total magnetization for each Landau level can be 
obtained by multiplying the above quantities by the Lan- 
dau level degeneracy L^/(27r/^), representing the total 
number of states with characteristic length Is that the 
surface area can accommodate. Since the magneti- 
zation contribution computed above for each eigenstate 
is an explicit function of its energy, we can perform the 
following integral to find the total magnetization density 
due to the electronic spins 

Ms = J de D{e)Ms{e)nF{e), (15) 

where M.s{s) = QsI^b^/ {^2^) is the magnetization of the 
eigenstate with energy e and np^s) = l/[e'^^~^)/^^^ + 1] 
is the Fermi-Dirac distribution function. The electronic 
density of states associated with the surface states is 
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FIG. 3: (Color online) Spin magnetization in units of (xo- 
Tesla) as a function of the magnetic field for (from left to 
right) Ao = -2, -1,0, 1,2 meV and ksT = O.lmeV. = 8 
(top panel) and 30 (bottom panel). 



D{e). For the Hamiltonian we used in the previous sec- 
tion it takes the following form 



D{e) 



1 



n>0,Q!=d= 



, (16) 



where Uc = (A^ — A^)/e^ is the Landau level index be- 
yond which the energy exceeds the cutoff energy A. The 
cutoff can be chosen to be the energy where the surface 
band becomes degenerate with the bulk bands and here 
we assume A = 300 meV. Using this density of states 
function to perform integration in Eq. ( 15 ) yields 



5o 



Xo 



-npisQ) + sgn(5o) 



E 

n=l 



-^F(gn) 



(17) 

where xo = {cQsIIb) /'^h. 

For a constant magnetic field Bq the magnetization 
of a single TI surface given by Eq. (17) shows an inter- 



esting behaviour as a function of chemical potential /i 
illustrated in Fig. |2j In order to avoid ambiguity asso- 
ciated with the high-energy cutoff we choose to display 
5Ms = Ms{T^ii) — M5(0,0), i.e. spin magnetization rel- 
ative to the neutrality point at T = 0. For the negative 




FIG. 4: (Color online) The color at each point (50, Ao) rep- 
resents the magnitude of the spin susceptibility in the linear 
response regime, i.e., x — Mg/Bo, in units of xo- The step 
(diagonal white line) is given by Eq. (19) and results from 
n = Landau lavel crossing the chemical potential. The 
discontinuity evident at Bo = reflects the fact that suscep- 
tibility X diverges as Bo ^ 0. We have assumed Qs = S and 
ksT = O.OlmeV in this graph. 



values of /i magnetization initially decreases reflecting the 
fact that the negative-energy surface states exhibit nega- 
tive spin polarization, as can be seen from Eq. (13). The 



large jump in SMs near fi = results from electrons fill- 
ing the fully spin-polarized n = Landau level. Further 
increase in ja results in increase of SMg now reflecting 
the fact that the positive-energy surface states exhibit 
positive spin polarization. 

The above behavior is unique to topological insulators 
for it results from the Landau level structure in a single 
Dirac point (or more generally an odd number thereof). 
In a TI the 'other' Dirac point is located on the opposite 
surface where the magnetic field points in the direction 
opposite relative to the surface normal, see Fig.[l] The 
contribution of this surface to SMs would be the same. 
We emphasize that for a relatively thick TI slab /3-NMR 
will be sensitive to a single surface facing the beam and 
the behaviour predicted here is in principle observable, 
except that continuous tuning of the chemical potential 
in a crystal might be difficult to achieve. 

A much more feasible experiment involves varying 
magnetic field Bq while keeping /i constant. We now show 
that a unique signature of the fully polarized n = Lan- 
dau level still exists in samples with intrinsic magnetic 
ordering, when the chemical potential resides in the gap 
(i.e. both the bulk and the surface are insulating in the 
absence of the field). In his situation we can set /i = in 
our model. Now consider the effect of the applied mag- 
netic field. Since the gap between the adjacent energy 
levels in which the Fermi energy is located is fairly large 
(~ 100 K for a IT field), at sufficiently low temperatures 
we can replace the Fermi function in Eq. (17) by the step 
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FIG. 5: (Color online) (a) Schematic behaviour of Ao, i.e., the 
intrinsic magnetization gap on the surface in a system with 
a slab geometry shown in Fig. [l] versus the applied magnetic 
field, (b) The energy of the n = Landau level versus the 
magnetic field. Starting at the point 1 and by decreasing the 
magnetic field gradually, the energy would follow the path 
shown by the blue curve, i.e., 1 ^ 2 ^ 3 ^ 4, this happens 
if Ao is described by the blue part of the cycle in (a). Now 
by increasing the magnetic field from a negative value, cor- 
responding to the point 4, the energy would follow the path 
given by the red curve, i.e., 4 ^ 2 ^ 3 ^ 1, since this time 
Ao would be given by the red curve in the hysteresis cycle. 



function riF^e) Q{—e) and write 



Bo 



Xo 



-6(-£o)-sgn(5o)£ 



(18) 



We have assumed that the chemical potential remains 
pinned at zero energy and does not change as we tune 
the magnetic field. With these assumptions the magne- 
tization has an interesting discontinuity at a finite mag- 
netic field. The discontinuity in Ms occurs since tuning 
the magnetic field forces the fully spin polarized Landau 
level energy to evolve according to Eq. Q and cross the 
chemical potential. The critical magnetic field at which 
the jump happens is given by 



Br 



-Ao. 



(19) 



Assuming that the intrinsic gap Aq remains independent 
of Bq we plot the resulting magnetization in Fig. |4j This 
predicted behaviour could be employed to experimentally 
detect the intrinsic magnetization gap in the surface of 
a magnetically doped TI and measure the size of Aq 
through Eq. (19). Although this signature only occurs 



when the chemical potential lies very close to the sur- 
face state Dirac point we note that magnetically doped 
samples satisfying this requirement have been grown and 
studied. 

It is important to note that Aq is in general not inde- 
pendent of Bq since magnetic moments of dopants will 
tend to align with the applied field. Thus, like in a fer- 
romagnet, there will be a hysteresis effect whose features 
will depend on the material details^^. One possible sce- 
nario is shown schematically in Fig. |5] In the case when 
Ao depends on the field the equation ([l9| continues to 
hold but must now be viewed as an implicit equation for 
the critical field Be. 



Another question that arises has to do with the origin 
of the electrons that fill the n = Landau level upon 
changing the field through Be. One may wonder where 
the extra electrons come from in a fully gapped isolated 
system. The answer lies in the side surfaces which under 
generic conditions remain gapless and act as a reservoir 
of electrons. The model we consider here does not cap- 
ture these states but they reflect themselves in solutions 
of the Hamiltonian which are not normalizable in an infi- 
nite system. These are fully spin polarized with opposite 
energy and spin direction. In fact they are the particle- 
hole conjugates of the fully spin polarized LL given in Eq. 
([8|. Taking into account these electronic states and the 
fact that it is more favourable for electrons with higher 
energy to be transferred to the negative energy states the 
counting problem can be resolved. 



IV. THE KNIGHT SHIFT 

We now outline how a /3-NMR experiment can in prin- 
ciple be used to probe some of the physics discussed in 
the previous section. The valence and conduction elec- 
trons in a metal posses magnetic moments arising from 
both their orbital motion and their spin degrees of free- 
dom. Nuclear magnetic resonance technique can be used 
as a probe of the spin part of the total magnetization 
in the presence of the magnetic field by measuring the 
relative shift in the nuclear resonance peak with respect 
to the same resonance peak in a reference insulating sys- 
tem. This effect, which is due to the interaction between 
electronic spins and those of the nuclei, is known in the 
literature as the Knight shift and has been extensively 
studied in both theory and experiment.!^ 

The mobile electrons in a metal interact with the 
nuclei in their proximity and the Knight shift in the 
resonance peak of these nuclei can be described by a 
local Fermi contact interaction term given by i^int = 
— ^IIq- /J^iS{r — Ri) where Ri is the position of the 
iih. nucleus and iii = jN^i is its magnetic moment. The 
magnitude of 7Ar, the gyromagnetic ratio, depends on 
the nucleus quantum state. The nucleus total spin, 
couples to the applied magnetic field and therefore the 
position of the peak depends on the magnitude of the 
total magnetic field experienced by the nucleus which 
has a contribution due to the interaction with electrons. 
It turns out that this shift is proportional to the spin 
susceptibility. The constant of proportionality, known 
as the hyperfine coupling, can be computed using first 
principle calculations for the implanted nuclei. On the 
other hand, if we assume that the presence of the nuclei 
does not significantly alter the electronic states, then it is 
possible to approximate the shift for them by taking the 
expectation value of the aforementioned interaction term 
using the unperturbed electronic states. This is the low- 
est order approximation in the perturbative treatment 
of the interaction term. For metallic systems with spin 
degenerate bands this shift is proportional to the spin 



6 



susceptibility as it can been seen from a simple calcula- 
tion considering the fact that the spatial and spin degrees 
of freedom are uncorrelatecP^. 

The spin-momentum locking on the surface of a TI 
along with the energy dependence of the penetration 
depth can in principle change the above simple physics. 
Since it is not possible anymore to separate spin and 
orbital degrees of freedom, one might question the valid- 
ity of the linear relation between the Knight shift and 
the spin susceptibility. We devote the rest of this sec- 
tion to addressing this issue by considering a very sim- 
ple model. We assume that the nuclei do not alter the 
electronic states around them. It is important to note 
that this assumption may break down for the implanted 
nuclei if they modify the electronic states around them 
significantly and computing the Knight shift would then 
require a first principle calculations. 

The field experienced by the ith nucleus due to the 
interaction with the proximate electrons is given by 



SB,= 



Sir 



-fe{o-S{r - Ri))T, 



(20) 



where (...)t is the expectation value of over electronic 
states at temperature T. Therefore, the effective Hamil- 
tonian for the ensemble of nuclei takes the form 



iff =-?i7iv^/,-(Bo + 5Bi). 



(21) 



This way, zth nucleus would have a resonance peak 
= In {Bo + 5 Biz)' The Knight shift is then defined by 
comparing the resonance frequency with the frequency in 
a similar material without these electronic states 



Ki = 



UJQ 



Bo + SBi^ 



Bo 



Bo 



SB,, 
Bo 



(22) 



The shift in the resonance peak of the nuclei ensemble 
is the average of the knight shift from each individual 
nucleus and is given by 



K 



1 ^SB, 



N 



Bo 



(23) 



where TV is the number of the implanted nuclei. Using 
the electronic eigenstates given in Eq. (|6][8| we get the 
following expression for K 



3NBo 



E E ic(^oi'(ic,ti'(^^o-ic,;i'(i^±.)). 



(24) 

Here |V^^(z)p appears as a factor in the realistic 3D elec- 
tronic wave functions of the surface electrons reflecting 
the fact that the electrons have an energy dependent pen- 
etration depth into the bulk. The nuclei implanted in the 
system have a spatial probablity distribution Pnuc(^5 ''"x), 
which depends on the energy and diameter of the beam 
of the ions used in the /3-NMR experiment. If the dis- 
tribution function is known, we can replace the above 



summation with a 3D integral over the crystal volume 
^E ^ j dzd^r^P^^,{z,r^). (25) 

i 

Assuming that the distribution is uniform in the plane of 
the surface, i.e., P^uc{z^r±_) = P(^), we get 



K 



Stt 



(26) 



where we have performed the integration over the in- 
plane degrees of freedom and replaced the summation 
over kx with the LL degeneracy. We have also defined 
the nth LL weight, as 



(27) 



Now if we assume that different LL have the same pen- 
etration depth we have = fo for all n < Uc and we 
obtain 



K 



Snfo 1 
3^0 1% 



occ 

E- 

n,Q! 



Stt/o 



(28) 



We thus recover the linear proportionality of the 
Knight shift to the surface electronic spin susceptibility 
under reasonable assumptions. Note that if we relax the 
assumption that different LLs can now have different 
penetration depths, then the Knight shift would no 
longer be linearly proportional to the total spin suscep- 
tibility. Instead, it would be a weighted superposition 
of contributions from each individual LL to the spin 
susceptibility. Nevertheless, the Knight shift will still 
display the interesting behavior discussed in this study 
as long as is a reasonably slowly varying function of 
n. We should emphasize once again that although above 
considerations elaborate on the differences caused by 
the spin-momentum locking and the energy dependent 
penetration depth, they do not take into account the 
fact that the electronic wave-functions could be altered 
by the presence of the implanted nuclei and the hy- 
bridization with the adjacent nuclei. 



V. CONCLUSIONS 

The magnetic response of spins of the Dirac-like elec- 
trons on the surface of a topological insulator shows inter- 
esting features both in the absence and the presence of an 
intrinsic gap Aq. When the surface states are gapped ow- 
ing to the time-reversal breaking perturbation (i.e. due 
to magnetic doping) the n = Landau level, which is 
fully spin polarized, can have positive or negative energy 
depending on the sign of the intrinsic gap Aq relative to 
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that of the apphed magnetic field. It will therefore be 
completely filled or empty when the chemical potential 
is tuned to zero energy. Our study shows that this struc- 
ture results in an observable jump in the spin susceptibil- 
ity, mesurable e.g. through the /3-NMR Knight shift, as 
one tunes the applied magnetic field through the critical 
value Be given in Eq. (19). The effect may be used as a 



means to measure the magnitude of the intrinsic gap on 
the surface of a magnetically doped topological insulator 
if /3-NMR or another surface-sensitive technique could 
capture the magnetic response of the electronic spins on 



the surface. This behaviour is a unique feature of the 
topological insulator exotic surface states closely related 
to the special form of the spin- momentum entanglement. 
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